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O . Abstract 
(N 

Q^' The M0ller tetrad gravitational energy-momentum expression was recently 

D . evaluated for a small vacuum region using orthonormal frames adapted to Rie- 

mann normal coordinates. However the result was not proportional to the 
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' Bel- Robinson tensor Bapfj,u- Treating a modified quasilocal expressions in a 

^ ■ similar way, we found one unique combination that gives a multiple of B^p^u 
which provides a non-negative gravitational energy- momentum in the small 

qh| sphere approximation. Moreover, in addition to B^js^u^ we found a certain 

^ ■ tensor Sap^u + Ka/s^u which gives the same "energy-momentum" density in 

bX), vacuum. Using this tensor combination, we obtained an infinite set of solutions 

I that provides a positive gravitational energy within the same limit. 

> : 

^ ; 1 Introduction 

00 ■ 

^ ! Finding an appropriate quasilocal expression for the gravitational energy is an open 

^ I problem in general relativity. How do we even know that gravitational energy exists? 

I Tidal friction from the tidal force of the Moon on the Earth leads to the slowing of 

O ' the rate of Earth rotation. The lengthing of the day is the indirect evidence for the 

'i> ■ gravitational energy. 

^ , All matter and all other interaction fields, through their energy-momentum den- 

^ I sity, act as the source of gravity. These sources exchange energy-momentum with the 

■ - - ' gravitational field; all attempts to identify an energy-momentum density for gravity 

itself led to reference frame dependent quantities (i.e. pseudotensors) , a reflection of 
the fact that energy-momentum of an isolated gravitating system is inherently non- 
local. This feature can be understood in terms of the equivalence principle: gravity 
cannot be detected at a point. However, one can get around this difficulty using 
the idea of quasilocal energy-momentum, i.e., associated with a closed 2-surface sur- 
rounding a region [1] . 

The Bel-Robinson tensor has a desirable property as it provides a non-negative 
value, although it cannot be interpreted as a "stress energy" for gravity directly 
because it has the wrong dimension. The energy density has dimension cm~^, while 
Bap^yf^t^t^t^ has dimension cm~^, where is the timelike unit normal. However, 
the small sphere region limit can resolve this mismatch. In the small sphere limit, 
the quasilocal expression for the energy-momentum density should be a multiple of 
r^B^QQQ ~ ^'nr^ {j-"^ B ^qqq) [1], where r is the radius of the Euclidean volume 3-ball. 
Indeed r^5oooo matches the energy density dimension. 

A positivity of gravitational energy proof was obtained in orthonormal frames 
[2]. Success on a large scale automatically implies the positivity on the small region 



limit. Recently M0ller's tetrad gravitational energy-momentum expression [3] was 
evaluated for a small vacuum region using orthonormal frames adapted to Riemann 
normal coordinates. This result for the gravitational energy in the small sphere 
approximation is not positive definite. Treating a modified quasilocal expression [1] 
in a similar way, we found one unique combination that gives a multiple of the Bel- 
Robinson tensor which means that the gravitational energy is definitely non-negative. 

The Bel-Robinson tensor component -Boooo or -Boo/' gives the non-negative grav- 
itational "energy" density. We discovered that the sum of the tensor components 
5*0000 + Kqoqo or Soo/ + Kqqi offers the same "energy" density value. Based on this 
criteria, we also obtained an infinite set of solutions that provide a positive gravi- 
tational energy within the same small sphere limit using certain modified quasilocal 
expressions. 



2 Orthonormal frames and quadratic curvature 

The quasilocal quantities for small regions can be studied by Taylor expanding the 
Hamiltonian, including the divergence of its boundary term in a small spatial region 
surrounding a point. The reference is the fiat space geometry at this origin. The 
orthonormal frame satisfies 

e^(o) = 9^6^(0) = 0, rv(o) = 0, (1) 

2 a ^ a a a ^ a 

dijC^'aiO) = —-{R"iaj + R^jai), djT"i3i{0) = -R"i3ji, (2) 

where the Latin letters refer to coordinate frames (holonomic frames) and Greek letter 
means the orthonormal frames (non-holonomic frames). 
The Bel-Robinson [5] tensor is defined as 

Ba/Sfiu '■= RaXnaRp^v" + *RaXtia * R/3^v'^ 

RaXpcr Rp^u" + RcxXuaRp^ — -QapRxap^.R'^'''' (3) 

where the dual curvature is ^Rappu = ^^a/sxaR^'^ fiu- Furthermore, the tensors Sai3p.u 
[5] and Kap^u in vacuum are defined as 

Sa/3p.u '■= Rap.XaRpu^'^ + RauXaRpp''^^ + -^Qa/Sg puRxaprR^'^'^ , (4) 

3 

Kal3pv '■= RaXPaR^i^u' + RcxXflaRp^ p^ — -Qapd pvRxaprR'^'^''^ ■ (5) 

O 

The identity in vacuum RxappR'^^^u = ^dpuRxaprR^'^^^ is useful. 

It is known that the Bel-Robinson tensor is completely symmetric, we have found 
the identity 

= Sa/3pu + Sap/3u + S^uPp + K^Ppu + Ki^^p^ + K^uPp- (6) 



Technically, the alternative form can be written as 

The tensors Sap^i, and Kapnu are both symmetric at the first pair and last pair of 
indices. 

It turns out that in vacuum the small region energy density has a RNC Taylor 
series expansion of the form 

t 1/ — t i/ijX X , (^) 

the corresponding energy-momentum is 

^ fiab I X X dj X 



t=0 
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— ^ ixa -j^g ) (9) 

where a, 6 = 1, 2, 3. Note that = - t°^o°- In particular = ^Voo as 

Bapnv is completely traceless. More covariantly 

^O/iOO '■— tanlS-yt^t^f^ ■ (10) 

The "energy-momentum" associated with the Bel-Robinson tensor is 

where the electric part Eab and magnetic part Hab are defined in terms of the Weyl 
tensor as follows: 

Eab '■— Cambnt^t^j H^b '■— *Cambnt^'t^- (12) 

Moreover for Sa/Sfiv + KaPfn^, we have the following identity related with the Bel- 
Robinson tensor components 

'S'^ooo + Kfiooo = B^ooo = S^o/ = -S'^^^o/ + -f^/io/- (13) 

It means that <S'^ooo + -f^iuooo or S^oi'' + Khq/ have the same physical quantities as B^ooq 
or Bi^oiK 

3 Modified quasilocal boundary expressions 

For a first order Lagrangian density: 

C = dqAp- A{q,p) (14) 



where q is an f-form, p is a 3-form and A is the potential. The modified quasilocal 
expressions |H E| can be briefiy summarized as follows 

Bc^^c^iN) = Bp{N) + Ci^ArAg A + ec2Ag A in^P, (15) 

where Bp{N) = i^q A Ap — eAq A znP, e = (—1)-^ with /-form, ci and C2 are real 
numbers, Aq = q — q, Ap = p — p, q and p are the background reference values. For 
GR 

£ = i?"^Ar//, (16) 

so let 

q^T%, p^^Va". (17) 

Allowing for the background connection T^'p = 0, then (ITS!) becomes 

2/€ i3,„,,(Ar) = r"^ A iNVa^ + ci^^r°^ A A7]J - C2r"^ A InAtjJ. (18) 

When (ci,C2) = (0,0), (0,1), (1,0) and (1,1), the four quasilocal boundary expres- 
sions with the simplest boundary conditions are 

Bp(0,0) = OAz^v///, (19) 

Sd(0,l) = OA^Tvr//, (20) 

i3e(l,0) = r%A2^r// + ^^r-^A Ar^/, (21) 

-B,(l,l) = r% A z^rj/ + ^jvF"^ A Ar//. (22) 



In terms of the superpotential, rewrite the modified quasilocal expressions (1181) as 

2k B,,,,,{N) = ~]^N^ + cit/^[^^^] - C2f/;.[*^1} e,,-. (23) 

The tetrad teleparallel gauge current expression in orthonormal frames is 

UM = -egP^T%m5ZV'^,e\e%, (24) 



and in RNC 



cif/^[*^l = -^e/'^i?"^^^i?VAeVe\xVx''5^;^ + 0(x^), (25) 
C2t/^[*^'] = ^eg^"R^p,xR''i.re\eKx^x^x^5Z% + 0{x^). (26) 

It should be noted that both the tetrad teleparallel gauge current U^'^^^ and the 
associated energy-momentum density dj{U^^^^^) is a tensor. In contrast, the M0ller 
1961 expression mUh^^^^ is a tensor but the corresponding energy-momentum density 
<9j(Aff^/i '*■'') is not a tensor. Precisely it is a pseudotensor, which means it depends on 
the coordinates in a non-covariant way. As mUh^'^^ = e^'hU^^'^l modify the superpo- 
tential of ([23]) to 

2«;ZY;,fo] = A^f/,[^^'] + e\ (cif/^[*^] - C2f/^I^^^1) . (27) 



Then the corresponding pseudotensor becomes locally in a small region in RNC 

= e2GH' 

+^ {2(2 - 3c2)Bh'^. - (1 - 3ci + 3C2)5/,V + 2(ci - 2c2)ir„V} 

+0(Ricci,a;) +C(x^). (28) 

Regarding whether the above expression is good for the gravitational energy, there 
are three limits we can consider. They are inside matter (interior mass density), at 
spatial infinity (ADM mass energy) and in vacuum (positive gravitational energy). 
The first two tests are relatively mild in general, but the last one is not, because it is 
very sensitive if we insist to obtain the Bel- Robinson tensor. For example, the M0ller 
1961 expression fulfills the first two tests while it fails the third examination [3]. 
Test (i): Inside matter. The energy density inside matter at the origin is 

S = -to°(0) = = -To°(0) = p. (29) 

K 

Test (ii): At spatial infinity. The total energy of the pseudotensor agrees with the 
ADM mass formula |7j 

E = yJ = ^ Jim Y: / {h,,, - h\,) N^dA, (30) 

where the integrals are taken over a sphere of constant r and A^-' = x-' /r is the outward 
normal to this sphere. 

Test (iii): In vacuum. Consider (1251) by eliminating the tensors Sa/s^u and K^/s^u 
when (ci, C2) = (|, |), the gravitational energy-momentum density in the small sphere 
region limit is 

ta^ = ———BJ^Ci^X^X . (31) 

12k 

This is the first desired result we found in vacuum. It is proportional to the Bel- 
Robinson tensor which is an invariant strength measurement of the non-negative 
gravitational energy density within a very small region. Explicitly 

-Boooo = EabE""^ + HabH'^^ > 0. (32) 

The second desired result for the non-negative gravitational energy-momentum in 
the small sphere limit is 

= --^ f {2(2 - 3c2)B^oij - (1 - 3ci + 3c2) + 2(ci - 2c2)K^oij} x'x^ 
2ci - 1 



X 
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provided ci > 1/2 and the unique combination ci + C2 = 1, which is the constraint 
such that the coefficients of 5*^0/ and K^q/ are the same. There is an infinite set of 
solutions because of the constant Ci. When (ci, C2) = (1, 0), the quasilocal expression 
for this set is Bc{l,0) as mentioned in (^1]) which has a simple boundary condition. 



4 Conclusion 



Once the positivity energy proof is achieved for some particular expression, the small 
sphere limit is guaranteed for the positive gravitational energy calculation. Recently, 
the M0ller tetrad gravitational energy-momentum expression was evaluated for a 
small vacuum region using orthonormal frames adapted to Riemann normal coordi- 
nates [3]. Treating a modified quasilocal expressions in a similar way, we found one 
unique combination that gives a multiple of the Bel-Robinson tensor. 

The components of the Bel-Robinson tensor -Bqooo = -Boo/ gives the non-negative 
"energy" density which has a nice property for the gravitational field. However, 
besides this tensor, we found that 5*0000 + -f^oooo = Sqqi^ + Kqqi'- gives the same physical 
value. Based on this property, we also obtained an infinite set of solutions that provide 
the positive gravitational energy within the same region limit using certain modified 
quasilocal expressions. 
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